Transition to quantum turbulence in a finite size superfluid 
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A novel concept of quantum turbulence in finite size superfluids, such as trapped bosonic atoms, 
is discussed. We have used an atomic 87 Rb BEC to study the emergence of this phenomenon. In 
our experiment, the transition to the quantum turbulent regime is characterized by a tangled vortex 
lines formation, controlled by the amplitude and time duration of the excitation produced by an 
external oscillating field. A simple model is suggested to account for the experimental observations. 
The transition from the non-turbulent to the turbulent regime is a rather gradual crossover. But it 
takes place in a sharp enough way, allowing for the definition of an effective critical line separating 
the regimes. Quantum turbulence emerging in a finite-size superfluid may be a new idea helpful for 
revealing important features associated to turbulence, a more general and broad phenomenon. 
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I. INTRODUCTION 



Quantum Turbulence (QT) is a phenomenon related 
to the vortex dynamics in superfluids and it can be real- 
ized in many different ways. In the case of liquid helium, 
below the A-point, moving grids and vibrating objects 
can generate a fully tangled configuration of quantized 
vortices which characterizes QT [J Q ■ Within the con- 
text of low temperature physics, QT has been studied 
since its discovery over fifty years ago Q. One of the 
main motivations of studying QT is to establish its re- 
lation with turbulence in classical fluids, where there is 
no requirement for the vortex angular momentum quan- 
tization. For many years, QT could only be studied in 
liquid helium ( 4 He and 3 He) However, recently [f| 
Bose-Einstein Condensate (BEC) of trapped gases has 
provided the main ingredients to study QT in a more 
simple superfluids. The vortex nucleation in a BEC can 
be produced by introducing a rotation in the trapped 
cloud 6-8}. When quantized vortices are generated, they 
arrange to form the Abrikosov lattices Q- These crys- 
talline structures of quantized vortices result from mutual 
interactions of the rotating field producing a repulsion 
between vortex pairs which is balanced by the trapping 
potential. In such structures, with a collection of vortices 
having the same direction and circulation sign, QT can- 
not spontaneously occur. Since the main characteristic of 
QT is a spatial tangled distribution of quantized vortices, 
one must have vortex lines distributed in many spatial 
directions to reach such a configuration. One method 
to achieve this was proposed by Kobayashi and Tsubota 
10]. In their proposal, combined rotations around two 
orthogonal axes induce the nucleation of quantized vor- 
tex lines in orthogonal directions with a clear evolution 
to QT. In a recent work [ll[ a variant of the procedure 
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suggested in Ref. [10| was implemented. Henn et al. [Tl( 
demonstrated that a special type of oscillatory excitation 
imposed on a BEC generates vortices. When vortices 
and anti- vortices are formed and proliferate through the 
sample, the emergence of turbulence is observed [f| as 
a configuration of tangled vortices. In fact, such oscil- 
latory excitations generate vortex-antivortex pairs [l2j]. 
The method of creating vortices by means of oscillating 
external fields is a particular case of the general method 
of creating coherent topological modes by such oscillating 
fields [iMl]. 

Quantized vortices inside a condensed atomic cloud 
are the necessary ingredients to produce quantum tur- 
bulence. Once the cloud is filled with vortices and anti- 
vortices distributed in various directions, but not ar- 
ranged in a lattice, turbulence should naturally be es- 
tablished. The production of vortices in a BEC is there- 
fore the first step. The majority of experimental groups 
have produced separate vortices and vortex arrays by in- 
troducing a single- axis rotation in a BEC. The stirring 
technique has been used with great success. In these ex- 
periments a similarity with experiments performed with 
4 He [HI could be directly observed. For large size BECs 
produced by the MIT group Q a large number of vortices 
with the same circulation could be produced and config- 
urations like Abrikosov lattices (originally observed for 
the type-II superconductors) were clearly seen. As much 
as 130 vortices could be accommodated in the large MIT 
BEC. In such experiments, however, no tangled configu- 
rations could be observed, but only crystalline structures 
were produced. 

In this communication, we analyze the process of gen- 
erating QT in BEC, focusing on the consequences of 
having a finite-size sample which constitutes an impor- 
tant property of superfluids originating from cold trapped 
atoms. In this case, the formed vortices spread inside the 
trap until turbulence is developed. We analyze this ef- 
fect in terms of the amplitude and time of excitation by 
the oscillatory perturbation. We suggest a simple model 



demonstrating the existence of an effective transition line 
between non-turbulent and turbulent regimes. A com- 
parison with experimental data is presented, showing a 
good qualitative agreement. We start with a brief de- 
scription of the emergence of QT in a BEC, followed by 
a general analysis and a comparison with experimental 
observations. Finally, arguments are presented concern- 
ing the importance of the phenomenon of quantum tur- 
bulence in finite superfluids. 



II. EXPERIMENTAL OBSERVATION 

The vortex nucleation can be done by imposing an ex- 
ternal oscillating field on the condensed trapped atomic 
cloud, producing an effect equivalent to rotation. In- 
deed, those oscillations generate vortex-antivortex pairs 
|13l - [l5| . Vortex nucleation is expected to take place when 
the perturbation provides enough energy. 

Our system is composed of 87 Rb atoms forming a 
BEC described earlier p| [ll|. We used a combina- 
tion of oscillations introduced by the superposition of 
coils to the conventional trapping field. The anisotropic 
harmonic potential of a quadrupole-Ioffe-configuration 
(QUIC) trap is approximately given by V = ^m(uipP 2 + 
ui 2 z 2 ). In our case, oj p /oj z — 9, which corresponds to 
a cigar shaped trap. Superimposing off-axis coils on this 
trap modifies the potential so that the long axis is rotated 
while the minimum is displaced. This type of excitation 
can produce different numbers of vortices, depending on 
the oscillation amplitude and excitation time This 
type of excitation should create both vortices and anti- 
vortices. Observations of three vortices inside an atomic 
cloud provide good evidence for the existence of vortex 
and anti- vortex pairs, as discussed elsewhere (T^ |. Since 
the vortex and anti-vortex pairs are not exactly parallel 
to each other (oscillations are not confined to a plane) 
the vortex and anti-vortex are expected to live longer 
than if they would be parallel. We have fixed the excita- 
tion frequency of 200 Hz and considered a range of am- 
plitude and time intervals of the applied external field. 
The results are observed after a time-of-flight (TOF) of 
20 ms. The overall cloud characteristics are determined 
by the trapping harmonic potential, whose frequencies 
are f z = 23 Hz and f p = 207 Hz. The cigar shaped BEC 
contains about 2 x I0 5 atoms with a typical condensate 
fraction ranging from 40 to 70%. 

As the oscillation amplitude/time increases, the vor- 
tices start to be nucleated in agreement with Ref. fllj j . 
The number of vortices in the cloud is counted as the 
clearly distinguishable dark regions in the absorption im- 
age after TOF and it changes with the amplitude and 
time interval of the excitation [17j. We have varied the 
amplitude from to 200 mG /cm and the excitation time 
up to 60 ms. We are interested in studying the crossover 
region from the observation of vortices to the turbu- 
lent cloud. The region in the diagram, separating both 
regimes, as well as the typical images characterizing the 
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FIG. 1: Atomic optical density images of: (a) non-turbulent 
cloud with well-defined separate vortices and (b) turbulent 
cloud, where the partial absorption changes along the image 
due to the existence of tangled vortices. The images were 
taken after 15 ms of free expansion. 



so-called non-turbulent and turbulent regimes are shown 
in Figffland FigH 

In order to distinguish between the non-turbulent and 
turbulent regimes, we have adopted the following crite- 
rion. In FigQJa), we present a density profile, where one 
can clearly see the vortices as the dark regions (valleys) 
spread inside the cloud. In this case, the typical Thomas- 
Fermi aspect ratio inversion is observed during the TOF 
(free fall), as expected in a non-turbulent gas. On the 
other hand, the turbulent regime is characterized by tan- 
gled vortex lines spread all over the cloud. As a result, 
the dark region contrast essentially fades away, as is seen 
in the absorption images of Figfljb). Besides that, the 
behavior of a QT cloud during the TOF is different than 
that of a regular BEC. The axes aspect ratio is kept es- 
sentially unchanged The tangled vortices make the 
whole system isotropic [l8|]. Combining these character- 
istics, we define a turbulent cloud. This is the criterion for 
distinguishing between the turbulent and non-turbulent 
regime. 
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FIG. 2: Diagram, on the excitation amplitude-time plane, 
demonstrating the crossover between the non-turbulent 
atomic cloud, with well-defined separate vortices, and the tur- 
bulent cloud, with tangled vortices. The plotted line is based 
on Eq. (fjj considering halfway points with the same ampli- 
tude as experimental critical points. 



III. THEORETICAL MODEL 

According to Figf2J there exists a well-defined parame- 
ter region, where the non-turbulent regime evolves to the 
turbulent regime. To understand the observed behavior, 
we propose a simple model based on energy-balance ar- 
guments, which are quite general and do not depend on 
the actual mechanism type needed for the vortex cre- 
ation or its dynamics inside the cloud. The nucleation of 
vortices is due to the instability of collective excitations 
arising from the energy pumped by external oscillating 
perturbations [l8l [l9|. 

We start by noting that there should exist a certain 
energy amount that is necessary to pump into the super- 
fluid atomic cloud for the vortex formation. Following 
[20| , we write down the energy needed for the vortex nu- 
cleation as 
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where, £ = (y / 87ma s ) -1 is the healing length, n is the 
BEC peak density, a s is the s-wave scattering length, and 
lo is the vortex line length. We assume that the latter 
is approximately equal to the effective cloud's harmonic 
oscillator length, 
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Then, let R pU m P be the rate of the total energy pumped 
into the cloud by the external oscillating field, and r\ be 
the energy fraction converted to rotation. Therefore the 
total energy needed for the vortex formation can be writ- 
ten as 



TjRpump \t to 
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where t is the elapsed excitation time and to is the mini- 
mal time period needed for the first vortex creation. The 
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the time t of pumping, the number of vortices N vort is 
formed, then the energy balance implies that 
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and the expected number of vortices is 
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This gives us a good estimate for the number of vor- 
tices observed in the atomic cloud, as a function of the 
excitation elapsed time, which can be compared to the 
time dependence reported in [l7|. Here, we have not 
explicitly taken into account the vortex-antivortex anni- 
hilation, but this effect could be incorporated into the 
value of the coefficient r\. This estimate depends neither 
on the cloud motion during the excitation, nor on the 
presence of collective modes that certainly arise |2l| . As- 
suming that the turbulence onset takes place when the 
atomic cloud is heavily populated by vortices, and N vort t; 
is of the order of the characteristic trap size, we conclude 
that turbulence should arise when the number of vortices 



Near the frontier between the two regimes, the above 
equations result in the critical behavior of R pum pi gener- 
ating turbulence: 
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The energy pump rate, R pump , is proportional to the 
ratio between the external field amplitude, A, and the os- 
cillating (pump) frequency. Therefore, the critical pump- 
ing amplitude can be expressed as 



A c (t) = 



C 



t-tn 



(7) 



The borderline separating the non-turbulent and the 
turbulent regimes is shown in Figj2j For our experimental 
system, we have found C rs 1.6 ms(G/cm), to ~ 17ms, 
£ ps 0.06^m, lo ps 1.08/rai, and E vort ps 20nKxkB, 
as the characteristic values. With these quantities, we 
expect N vort ps 20 for the QT onset, which is in good 
agreement with the recent work [17j . where the turbu- 
lence is observed when the number of vortices is close to 
the value 20 found above. 



IV. ROLE OF THE FINITE SIZE OF ATOMIC 
SYSTEMS 

The existence of the critical borderline, separating the 
non-turbulent and turbulent regions in the diagram of 
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Fig[2] is closely related to the finite-size of the trapped 
superfluid system, whose effective size is given by the 
characteristic anisotropic harmonic trap length Iq. Quan- 
tum turbulence develops when the trapped atomic cloud 
becomes densely filled with vortices. Then the energy, 
pumped into the system, transforms not only into the 
newly formed vortices, but also to their rapid motion, 
with the formation of their tangled distribution, accom- 
panied by the appearance of reconnections and the for- 
mation of Kelvin waves. At this point, the absorption im- 
ages become hazy, which is a manifestation of the arising 
turbulence [22 1. 

The standard experiments with 4 He and 3 He do not 
fulfill the conditions, where finite-size effects would be- 
come crucially important. Hence, the superfluid clouds, 
formed by trapped atomic BECs, represent a whole new 
class of systems, whose properties essentially depend on 
finite-size effects, as well as on interactions. 

Effects, related to finite temperature, should also be 
important, though, at this time, we do not have enough 
data for performing the corresponding analysis. The 
pumped energy is, certainly, partially transformed into 
the thermal cloud, which is necessary to take into ac- 
count for a more detailed consideration [23[. It can also 
be that the losses, caused by the thermal cloud, could be 
responsible for the value of the delay time to. The exis- 
tence of the critical line in Fig J5] requires that t be larger 
than to- ln order to generate vortices at t ~ to, one needs 
a very large amplitude. But a strong or long pumping 
should produce a large admixture of the thermal fraction 

El. 



fluid formed by a Bose-Einstein condensate of trapped 
atoms. The transition region occurs on the amplitude- 
time parameter of the related excitation. The charac- 
ter of the observed transition region is closely connected 
with the fmitcness of the superfluid, since quantum tur- 
bulence arises when the sample becomes densely satu- 
rated with vortex lines. A simple model allows us to 
qualitatively understand the observed main features of 
the phenomenon. 

Further studies, both experimental and theoretical, 
are necessary for the better understanding of this phe- 
nomenon of quantum turbulence in finite-size superfluids. 
This concept of quantum turbulence in finite systems can 
also be applied to small 4 He droplets, though their exper- 
imental realization can be much more complicated than 
the creation of atomic clouds in traps. The possible exis- 
tence of finite-size effects in turbulent superfluids opens 
a novel direction in the investigation of finite systems, 
such as trapped atoms and liquid droplets. 

At low temperatures, fermions, depending on the sign 
of their interactions, can form either superfluid molecular 
BEC or paired superconductor-type fluid, both of which 
can exhibit superfluid properties @, It is, therefore, 
feasible to produce quantum turbulence in such systems 
and to observe finite-size effects in fcrmionic turbulent 
fluids, similar to those observed in trapped bosonic su- 
perfluids. 
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